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Abstract 

The lifetime of a disoriented chiral condensate formed within a heat bath of 
pions is calculated assuming temperatures and densities attainable at present 
and future heavy- ion colliders. A generalization of the reduction formula to 
include coherent states allows us to derive a formula for the decay rate. We 
predict the half-life to be between 4 and 7 fm/c, depending on the assumed 
pion density. We also calculate the lifetime in the presence of higher reso- 
nances and baryons, which shortens the lifetime by at most 20%. 

PACS numbers: 25.75.-q, ll.30.Rd, 12.38.Mh 
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Formation of hot and dense matter in heavy-ion colhsions has the possibihty of creating 
a phase where chiral symmetry is restored. As this matter cools and expands, the vacuum 
could relax into the "wrong" zero temperature ground state. Subsequent shifting of the 
vacuum back into alignment with the outside world could then lead to an excess of low mo- 
mentum pions [|I| in a single direction in isospin space 0. This excess is called a disoriented 
chiral condensate (DCC) and many studies have looked into whether its presence could be 
a signal of chiral symmetry restoration in heavy-ion collisions. 

Most of the discussion has centered around the details of formation of such regions . 
There has been a general consensus that these regions could be numerous and large enough 
with the conditions at RHIC and the LHC to be detected. While the idea of DCC-formation 
is appealing, the physics governing the chiral phase transition is not known well enough to 
make reliable predictions about the possible formation of these condensates. Experiment is 
needed to establish their existence. 

The ability to detect a DCC from hadronic observables in heavy-ion collisions depends 
on the condensate lifetime. An early estimate |^ gave a half-life of r ~ 3 fm/c. A more 
recent calculation has found a damping rate of 7 ~ 1 (fm/c)~^. This roughly corresponds 
to a half-life ofr~l/7~l fm/c, which would be short enough to jeopardize the definite- 
ness of the signal. This calculation was based on the 0(4) sigma model in the symmetric 
phase. However, after formation, the DCC lives in the phase of spontaneously broken chiral 
symmetry. Characteristic of this phase is the suppression of S-wave scattering among pions, 
which should protect the low momentum pion modes in the DCC. Indeed, an extension of 
Ref. to the broken phase produces smaller damping rates 0. 

The purpose of this letter is to provide a reliable estimate of the lifetime of a DCC state 
in the hadronic (chirally broken) phase. Defining the DCC to be a coherent state of pions 
with low momenta, we derive a general formula for the decay rate in the presence of other 
hadrons. The effect of higher resonances as well as that of baryons to the DCC lifetime is 
also estimated. Our calculation is constrained by data at every point possible. 

Assuming the formation of a DCC, interactions with the thermal heat bath can enhance 
or deplete the number of pions in the condensate. In a heavy-ion collision, pions are the 
most abundant thermal particles, and so their interactions with the DCC is expected to 
give the dominant contribution to the decay rate. The DCC can be written as a coherent 
state 0,| 

normalized to unity and for definiteness taken to consist of neutral pions, vr", created by 
a].. The Lorentz invariant measure dk = d^k/{2TT)^2Ek is weighted by the normalizable 
function rjk to give a momentum distribution which could evolve over time . The number 
of particles in the DCC is then given by = Jdk l^^fcp. 

We can get an estimate of the decay rate by considering two body interactions of the 
DCC pions with the surrounding thermal pions. Contributions come from vrvr scattering with 
either two thermal pions knocking one pion into the condensate or a DCC pion interacting 
with a thermal pion to escape from the condensate. There could also be individual scatterings 
with two pions knocked out or put into the condensate, but these should be suppressed due 
to the restricted phase space of the DCC. Using 
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{7]\rf\x)\r]) = Jdke-''^{r]\ak\7]) = jdke-'^^7]k, (2) 

a rederivation of the Lehmann-Symanzik-Zimmerman (LSZ) reduction formula taking into 
account the coherent state gives 
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The matrix element for tttt scattering (|7^^p) = |l]/|(2/ + l)|T^7rP is isospin averaged 
to account for only vr'^'s coming from the DCC and the factor of | accounts for identical 
particles in the final state after isospin considerations. The thermal weighting is given by 

Fu3 = /2/3(l + /l) - /l(l + /2)(1 + /s), (4) 

with fi = {exp{Ei/T) — 1)^^ representing the Bose-Einstein momentum distribution. The 
DCC pion energy is Eq = [{k2 + k^ — kiY + m^]^/^. 

The phase space correlations in Eq. (|^) can be simplified for DCC pions narrowly peaked 
near a single momentum k^. Then the square of the momentum distribution can be replaced 
by a momentum-conserving delta-function \rjk^^k^_k^\'^ A^(27r)'^ 2_E'o ^'^(^o + ki — k2 — fcs), 
giving the decay rate 

1 dN 



N dt 2Eo 



Jdk, dk2 dh Fu3 {\%n\') (27r) V(S k,) . (5) 



Although this formula pertains to two-particle interactions, it can easily be generalized to 
include more particles. To understand the various contributions to the decay rate, we first 
take the DCC pions to be at rest, ko = {jriT,, 0). We relax this condition in Eq. ([T3|) , allowing 
for a finite spread in the momenta in accordance with Eq. (^). 

This result was derived assuming the removal or addition of a pion to the condensate does 
not effect the DCC coherent state. This approximation is valid to better than 10% until the 
number of particles in the condensate decreases below ~ 5, at which point the existence 
of a macroscopic condensate is questionable anyway. Eq. @ can also be derived starting 
from the Boltzmann equation by assuming the DCC momentum distribution function obeys 
/o ~ 1 + /o, which again is appropriate for large enough and distributions peaked at zero 
momentum. 

We can now estimate the decay rate by replacing the tttt scattering matrix element by 
its experimental value in terms of the center of mass energy ^/s, three-momentum g, and 
scattering angle 6*, 

TL = 3277 E (2/ + l)PKcos 6) ^ e^'' sin . (6) 

Accounting for the three dominant phase shifts 6q, 6l, and 6q by using a parameterization 
of the data , the momentum integral in the decay rate can be readily evaluated. Defining 
the half-life r by 
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we find r = 8.9 fm/c for T = 150 MeV and r = 5.6 fm/c for T = 170 MeV. Considering the 
entire hadronic phase exists for about 10 — 20 fm/c, the DCC could live for a substantial 
fraction of the hadronic phase, allowing for a strong signal of its existence, e.g. in the dilepton 
channel 0. Even assuming three times the number of thermal pions, as predicted by some 
event generators [|T^], the half-life is still r = 3.8 fm/c at T = 150 MeV. 

It should be noted that the momentum dependence of the phase shifts is important to 
the result. Approximating the vrvr amplitude purely by threshold (g = 0) scattering lengths 
gives a half-life five times longer. Only if the range terms, which have a dependence, are 
also kept does the estimate come within 10% of the full phase shift result. The p channel 
(P-wave) contribution is merely a few percent effect, in agreement with estimates using a 
simple Breit-Wigner form. The main contribution to the decay rate comes from the S'-wave. 
This is because the other phase shifts are very small in the region where the Bose-Einstein 
distribution functions are appreciable. 

The p resonance could possibly play a role as a thermal particle though. We therefore 
look at vrp scattering. The dominant contribution comes from the formation of an ai(1240) 
meson in the s-channel, which we can model by a relativistic S'-wave Breit-Wigner matrix 
element 



q2 (,_^2^)2 + ,r^i(,) 



with center of mass three-momentum q and a momentum dependent width whose form is 
dictated by general decay of an axial-vector particle into a vector and pseudoscalar particle 



al ' 



q (3m^ + 
qai{5ml + qX 



-,,(.)= 400 ; . J. y^ MeV. (9) 



Adding this matrix element to Eq. gives a negligible effect, not changing the original tttt 
scattering half-life to the accuracy quoted. This can be seen in Fig. |l|, where the matrix 
elements for vrvr and vrp scattering are compared for various temperatures. 

The Boltzmann suppression of the ai pole along with threshold suppression from the 
width conspire to make Tip scattering a small effect. Other meson-meson processes not 
considered here, such as nai scattering, are therefore expected to be smaller. 

Baryons could also play a substantial role in reducing the lifetime of a DCC. Even at 
SPS energies (-Ecm = 20 GeV/nucleon) a large number of baryons are found in the central 
rapidity region. We therefore study the dependence of the half-life on nucleon density. The 
matrix element for vrA^ scattering is [|lT| 



{\T^^\') = 6An's E ^(2/ + 1) (\9iis, e)\' + Ms, e)\') (10) 
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with gi{s, 6) and hj{s, 9) denoting the spin non-flip and spin-flip amplitudes for each of the 

1 3 
2' 2' 



two possible isospins I = ^, |. These can be written in terms of phase shifts labeled by 



a = L2L2J 
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FIG. 1. The contribution to the decay rate, Eq. (||), in inverse fm/c from vrvr, tt/o, and ttN 
scattering keeping the ratio of pions to nucleons fixed at five. 

gi/2{s, 9) = Tsu + cos6' (2J>i3 + Tpii) + ... 

^i/2(s, 6) = jFpi3 — J^pii + . . . 

with similar expressions for the I = \ channel by direct substitution of the isospin indices. 
The ellipses refer to higher partial waves we will not consider here and 

JT^ = (gcot(5c,(g) -iqY^ • (11) 

Adding this matrix element to Eq. (|^), we must also replace the thermal weighting factor of 
this term by 

= hf^{l - fn - //"(I + /2)(1 - /s^) (12) 

with = (exp((i?j — fi)/T) + 1)^^. A chemical potential fi has been introduced to enforce 
a particular nucleon density. 

For a parameterization of the vrA^ phase shifts, we use the full relativistically improved 



A(1232) isobar model ||Tl]], which reproduces the experimental phase shifts very well up to 
the inelastic threshold. Adding this decay rate to those of vrvr and vrp scattering without a 
chemical potential gives only a few percent decrease in the half-life. But taking /i = 260 
MeV at T = 150 MeV to enforce a 5 : 1 ratio between pions and nucleons as observed at 



the SPS [|T2|, we end up with a half-life of r = 7.2 fm/c. The effect of the nucleons is about 
a 20% reduction in r. Still, this effect is small enough to imply experiments already taking 
place at the SPS should be capable of detecting remnants of DCC formation, provided of 
course that the conditions for DCC formation are met at these energies. 

A summary of our results is shown in Fig. |I[ At each temperature the chemical potential 
of the nucleons is chosen to enforce the 5 : 1 ratio between pions and nucleons. Since the effect 
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FIG. 2. The half-life r at T = 150 MeV for various pion to nucleon ratios. The dashed line is 
for three times the number of thermal pions and nucleons. 



of the nucleons is small, this result also applies for RHIC energies {Ecm = 200 GeV/nucleon) 
where the central rapidity pion to baryon ratio is probably even larger. 

We also show a plot of the dependence of half-life on nucleon density for T = 150 MeV 
in Fig. 1^. The ttN contribution becomes comparable with the tttt contribution when there 
is about one nucleon for every pion. However, we have so far assumed completely ther- 
malized pions and nucleons for this estimate, whereas some event generators indicate an 
enhancement of the occupation numbers by a factor of three [|l^. Taking this into account 
in our calculation gives the dashed line in Fig. 0, which shows the half-life is between 3 
and 3.5 fm/c for all the pion to nucleon ratios considered. If the over-population of particle 
spectra predicted by event generators is correct, it implies the DCC signal would be difficult 
to detect from hadronic observables, requiring an indirect means such as dileptons P]. 

Rather than taking the condensate to consist of only pions at rest, we can investigate the 
change in the lifetime for a finite spread in the DCC momenta. This will open up the phase 
space of the DCC and allow for further interactions with the heat bath. By introducing an 
integration over the Mandelstam variables s = {ko + ki)^ and t = (/cq — ^2)^ 
simplify Eq. (^) for ttP scattering (with P=7c,p, or A^) to 
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we can 
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FIG. 3. The average half-life (r) for T = 150 MeV as a function of the width of the DCC 
pion momentum distribution a from vrvr scattering alone (solid) and also including vrp and ttN 
scattering (dashed). 

= 3 = m]^, fn\ = m\ = m^, m\ = 2EqEz + u . 

with integration limits dictated by s > {m-,, + mpY, t < 0, Kf > 0, $ > 0, and Eq + Ei > 
E2 + mp. The thermal pion energy should also be restricted to be above the maximum DCC 
energy Eq. 

We then average the half-life obtained from Eq. ([l^) over the momentum distribution 
in the DCC as dictated by Eq. (|). Taking \rik\^ oc exp(— |kop/o"^), the average half-life is 
plotted as a function of the width a for T = 150 MeV in Fig. |^. The solid line is for the tttt 
scattering contribution alone and the dashed line also includes the np and irN contributions. 
The pion to nucleon ratio is kept fixed at 5 : 1. The decrease of the half-life from 8.9 fm/c 
to 7.2 fm/c for |ko| = was due to the nucleon contribution. As cr increases, P-wave 
7171 scattering and, to a lesser extent, 7ip scattering become more active, decreasing (r). 
Consequently, smaller DCC domains (characterized by a larger a) decay faster. 

In conclusion, we have derived a simple decay rate for a DCC and applied it assuming 
initial conditions at present and future heavy-ion colliders. In all stages of the calculation 
we have used input from data to constrain our results. The lifetime obtained is longer than 
previously estimated. 

We find interactions with the surrounding heat bath are dominated by vrvr scattering 
leading to r = 8.9 fm/c. Even with large nucleon densities, such as those that occur at the 
SPS, inclusion of vrA^ scattering reduces the half-life by at most 20%. Therefore models that 
only account for pion degrees of freedom should still give reasonable results. However, the 
factor of three estimated enhancement of thermal distributions assumed at the SPS is what 
most dramatically effects the half-life, reducing it to about 3-4 fm/c. This is probably not 
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long enough to observe DCCs in existing hadronic data. 

The decay of small DCC domains is accelerated by the contribution of P-wave tttt scat- 
tering as well as vrp scattering away from threshold, decreasing the average lifetime. With 
large domains formed in the absence of nucleons, as could occur at future colliders, the 
detection of DCCs looks more promising. 
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